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Abstract 

We study second order relativistic viscous hydrodynamics in 4-dimensional conformal field 
tiieories. We derive Kubo-type relations for second order hydrodynamic coefficients in terms of 
3-point stress tensor retarded correlators. For A/'=4 super Yang-Mills theory at strong coupling 
and at finite temperature we compute these stress tensor 3-point correlators, using AdS/CFT, 
by evaluating real-time cubic Witten diagrams in the AdS-Schwarzschild background. The 
small momentum expansion of the 3-point correlators in terms of first and second order hydro- 
dynamic coefficients is matched with the AdS result. We arrive at the same expressions for the 
hydrodynamic coefficients which multiply terms quadratic in the shear and vorticity tensors 
in the hydrodynamic expansion of the stress tensor as did Bhattacharyya, Hubeny, Minwalla 
and Rangamani [1]. Our method extends the results of Baier et al [2], and allows for a unified 
treatment of hydrodynamic coefficients, which are extracted from 2-, and now, 3-point retarded 
stress tensor correlators in the AdS-Schwarzschild background. 
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1 Introduction and Summary 



One of the most familiar applications of AdS/CFT [|3l, [|4] is computing linearized hydrodynamics 
for a variety of gauge theories with holographic duals (see [5J and references therein). More re- 
cently, second order hydrodynamic coefficients have been extracted using two different methods. 
Baier et al Q used Kubo-like formulae for 2-point stress tensor correlators to access a certain 
subset of second-order hydro coefficients of A/'=4 super Yang-Mills plasma. Their computation 
was made possible by a real-time finite temperature prescription for computing 2-point correlators 
[l6l 171. On the other hand, Bhattacharya et al [1] wrote a generalized black brane metric ansatz, 
dependent on the temperature and the black branes velocity viewed as collective fields, and turned 
the Einstein equations of motion solved perturbatively in the collective modes near the AdS bound- 
ary into equations of fluid dynamics. Specifically, the gravitational stress tensor, expanded near the 
boundary of AdS, took the form of a non-linear fluid dynamics stress tensor, with the various terms 
in the expansion in velocity and gradients multiplied by the corresponding hydrodynamic coeffi- 
cients. The subsets of hydro coefficients computed by [|2l and [|T|| agreed where they overlapped, 
and between them determined the full set of coefficients. By now, the second order hydrodynamic 
coefficients have been computed following [JJ in a variety of cases: at finite chemical potential 
dlllll, or in the presence of fundamental matter flO]. For recent review papers see filU and [|121 . 
The breakdown of second order hydrodynamics is investigated in fT3l. 

Our work is intended as a continuation of [|2l, where Kubo-like formulae are used in con- 
junction with higher-order stress tensor correlators to extract the second order hydro coefficients 
computed by jT]. We first derive the necessary Kubo relations using the method of Moore and 
Sohrabi [141^ then compute 3-point retarded stress tensor correlators in the hydrodynamic regime 
(i.e. in the limit of small momenta) in the AdS-Schwarzschild background which is dual to 
finite-temperature strongly coupled N=A super Yang-Mills in the limit of large number of col- 
ors Nc ^ 1. In other words, as opposed to [IJ, we do not deform the gravitational background, 
but instead compute higher-order correlators in the background. The problem of computing such 
higher-order real-time finite temperature correlators was solved in [15]. For momentum-space 
retarded 3-point stress tensor correlators this amounts to computing real-time Witten diagrams, 
depicted in Fig. lb, with three causal (two advanced and one retarded) graviton bulk-to-boundary 
propagators, joined at a bulk vertex which is integrated up to the black hole horizon (i.e. in the 
maximal causal diamondjl. 

To set up the problem, let us begin by revisiting some of the definitions and results of [|2l[T4l. 

The stress tensor of a conformal fluid can be written in terms of an equilibrium piece plus an 
extra term, which in the hydrodynamic regime can be series-expanded in gradients: 



'Readers should be warned that there is an error in the derivation of the Kubo relations for Ai, A2 and A3 in flT]. 
Specifically, they left out the (e + P)U^U^ + Pg^^ term in deriving their T^^. Also, = 0{h^) only in the static 
limit. 

^See also lfT6l for similar causaUty considerations and [ITTI for a different take on this subject. 
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(a) (b) 
Figure 1: Witten diagram for (a) 3-point correlator in imaginary time AdS-Schwarzschild and (b) 
retarded 3-point correlator of type raa with the boundary point x having the largest time; xi and 

X2 can have any time order. 



Here is the fluid's velocity, normalized such that 

U^U^g^, = -1, (3) 

and Yi^'' is defined such that U^Yi'^^ = 0. gfj_u is the background metric, whose curvature tensor is 
Rpfiu"^, and all derivatives are covariant with respect to the background metric, a^'^, il^'^ are the 
shear and vorticity tensors respectively: 

^ 2V^W^ = A^^A^"(Vpt/^ + VMp) - |A'^^A''"Vpf/, (4) 
= lA^^A^'^(Vpt/. - VMp) (5) 

where A^'' are transverse (to the fluid's velocity) projectors: 

Angular brackets denote transverse projection, followed by symmetrization and removal of a trace: 

Qif^u) ^ 1 A^PA^'^(Op, + O^p) - lA^'^A^'^Op^ (7) 

For a conformal fluid, energy density and pressure are related by the condition that the stress tensor 
is traceless, and so e = (rf — 1)P, where d is the number of space-time dimensions. 

In [14], the fluid's response to a small metric perturbation was computed by solving the stress 
tensor conservation law 

V^T^^ = 0, (8) 
together with the condition that the fluid describes a conformal theory, 

t; = 0, (9) 
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order-by-order in a double series expansion in the metric fluctuation and in gradients. This expan- 
sion of the stress tensor is then compared with 

+ ijd'xj d'yG^:t\<iz;x,y)h,^{x)K^{y) + ... (10) 

where Gra^.'.'a are retarded n-point correlators, with the measurement point z having the largest 
time. We assume the metric fluctuations hp^- vanish in the far past. To avoid clutter, for the 
rest of the paper we will suppress the ra, raa subscripts signifying retarded correlators^! with the 
understanding that all the correlators we compute are of this type. 

By identifying the two series expansions, namely the solution to (|8I9I) and the expansion in 
retarded correlators (ITOT i. one gains access to the hydrodynamic expansion coefficients in terms of 
causal stress tensor correlators. For example, the shear viscosity, i], can be computed from 2-point 
correlators expanded to linear order in gradients [19 1. Similarly, r]Tu and k can also be computed 
from 2-point correlators, expanded to second order in gradients |l2l. Ai, A2 and A3 require 3-point 
correlators, expanded to second order in gradients [JT4]. 

Since solving ([8]) is easier in momentum space, we will write the solution to ([8]) and ([91) as a 
series expansion in momenta. For a direct comparison with the response of the stress tensor we 
need then the Fourier transform of (HO): 



+ I J d% J d'q,j^^6\q - qi - q2)G^''\'"'\^Hq; -q,, -^2) V(9i)/^rc(92) + • • • 



(11) 



where we have used translation invariance to factor out momentum conservation delta- functions. 
In a 3-point function, there are two independent momenta, and the spatial momenta could point in 
different directions. However, for simplicity, we will assume that they both point in the z direction 
(i.e. we consider metric fluctuations which are independent of the x, y coordinates): 

qrf = 0, 0, fci), q^ = (0^2, 0, 0, k2). (12) 

The purpose of this paper is to obtain the remaining second order hydrodynamic coefficients Ai , A2 
and A3 from 3-point stress tensor correlators. We will find it simplest to study the sy-component 
{T^^)h and will extract second-order hydrodynamic coefficients from G"*^'' ' ': 

lim d^,d^^ lim G^'fl^-lf^ = -Ai + r^rn 
lim dk,d^, lim G^flJ/^lo- = -iA2 + ir/m 

lim dk,d,, lim G^^l°^l°^ = -iA3. (13) 

fcl->-0 oJi^O * 

fc2— >-0 aJ2->'0 



For a nice summary of the (r, a) notation see ifTsl . 
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We will find that the leading order AdS/CFT computation of the correlators leads to the same 
expressions for Ai, A2, A3 as those previously obtained by [ID: 

Ai = ^, A2 = , A3 = 0. (14) 

lb 8 



The paper is organized as follows. In Section 2 we derive the Kubo relations given in (fT3l) . 
More Kubo formulae can be found in Appendix |Bl In Section 3 we use AdS/CFT to evaluate the 
stress tensor correlators. We begin with a review of 2-point stress tensor correlators, and we com- 
ment how the hydrodynamic expansion of the retarded correlators and the AdS/CFT expressions 
are matched term-by-term in a gradient (small momenta) expansion. Then we present a similar re- 
sult for the 3-point retarded stress tensor correlators which appear in (fT3l) . with the final result for 
the second-order hydro coefficients given by (fT4] ). Technical details are relegated to appendices. 
For example, the r]Tu contribution to (fT3l) is derived in Appendix A. The graviton bulk-to-boundary 
propagators, expanded up to third order in momenta, are given in Appendix C. A few 2-point cor- 
relators are discussed in the text, but the rest of them are presented in Appendix D. Lastly, the 
on-shell (first, second and third order) gravitational action is given in Appendix E. 

We use the following conventions: five-dimensional tensors, which live in the AdS-Schwarzschild 
geometry, have indices given hy K, L, M, ■■■ = 0,1,2,3,5 = t, x,y, z,u; four-dimensional space- 
time indices are /i, z/, • • ■ = 0, 1, 2, 3; three-dimensional spatial indices are = 1, 2, 3 = 
X, y, z. 4-vectors are denoted by bold letters {x, q), and 3-vectors with an arrow (x). We use a bar 
to denote background values, in both the field theory and in the holographic dual. E.g. e is the field 
theory background energy density, and qmn denotes the background AdS-Schwarzschild metric. 
For simplicity of notation, for the rest of the paper we will also work in units where 2u:T = 1. 

2 Kubo formulae for Ai , A2 and A3 

2.1 The fluid velocity to leading order in the metric fluctuations and gradi- 
ents 

We begin by assuming that the fluid is initially in equilibrium in a flat space-time background, 
which is briefly distorted by a small gravitational perturbation: 

= r]^,u + h^^, e = e + 6e, P = P + SP (15) 

where h^,, {t, z) is a small metric fluctuation, and 5e, 5P are the induced energy density and pressure 
variations. In this section we will keep things general, and we will not enforce P = |e until the 
end. We take the fluid to be initially at rest 

U = U + 5U, f7^ = (l,0), (16) 

and so the background transverse projectors are simply 

^ij = ^ij^ ^oi = ^00 _ Q_ 

From the normalization condition f/^f/^ = — 1 we get the temporal component 

^7° = l + i/ioo + 0(/i'). (18) 
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The spatial components W are determined from the stress tensor conservation law: 

= V^T^" = [d^{e + P)]U^U'' + {d^P)g^'' + (e + P)V^(f/^f/'^) + V^n^^ (19) 

where 

V^,{UW) = d^{UW) + T^^^U^U"" + V^^U^U^ (20) 
To leading order in the metric fluctuations, this can be simplified using 

[d,{t + P)]UW + {d,P)g^' = d,{t + P)U^U'' + (9^P)r7^^ + 0{h^) (21) 

by noting that de and dP dSQO{h). 

First consider the u =±= y case and keep only terms through 0{h) in (fT9l) . Since tj^ = 
and d± = 0, the first two terms in (fT9] ) vanish at 0(h) by (|2T)) . Using (|20l) . and noting that both 

andrareO(/i), 

V^(f/^f/^) = 9^(t/''f/^) + T^^rU^U^ + 0{h^) 

= U^ + ho± + 0{h^), (22) 

where we use dots as shorthand for time derivatives. Now use the fact that 

t/± = g±uU'' = g^oU' + g^JJ^ + 0(/i2) = h^^ + + 0{h') (23) 

to rewrite this as 

V^iWU^) = U± + 0{h^). (24) 

Putting it all together (along with the fact that li^" is 0{h) and so VII ~ 911), the transverse 
case of the conservation law (fT9l ) is 

{t + P)U^ + d^U^^ = 0{h^). (25) 

This simplifies to 

{l+P)U^ + d,W^ = 0{e), (26) 

since H.^^ has only spatial components at and d± vanishes. 

Our goal will be to work to leading order in frequencies u and leading order in spatial momenta 
k, without assuming anything about the relative size of u and |A;|. (That is, we want to be able to 
handle w ~ A;^ and w ~ |A;| and a; = and A; = on an equal footing.) In (|26l) . we have explicitly 
kept the term which is leading order in time derivatives of U . We now need the term that is leading 
order in spatial derivatives of U . That comes from the first-order hydro terms in dilV-^. 

So consider the first order hydro expansion of Yi^^": 

= -r7A^-(V,f/;3 + V^f/JA^'^ + (fr^ - C)A^'^A"^V„t/;3, (27) 
where ( is the bulk viscosity and ^ = for a conformal theory. At first order in h. 

Kit = -vi^iUj + V,Ui) + ilf] - C)6,,VkUk + 0{h^), (28) 
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and in particular 

n^-^ = -vidzUi. + doK^ - dM + 0(/i2). (29) 

Putting this together with (|26|) . we get 

(e + P)U^- vK - vdz{k± - K^) + d,{2nd order hydro effects) = 0{h^), (30) 
where primes are shorthand for z-derivatives. The solution is 

Dr,k{kho± + ujh,±) ^ 
= ^ + OiuU^, kU^, M, (31) 

where 

is the diffusion constant associated with shear viscosity. Here (co, 0, 0, k) is the four-momentum 
associated with the factors of h. The 0{ujU_i_, kUj_) corrections listed in (1311 represent corrections 
due to second-order hydro, which are suppressed by additional factors of momenta^. Note that U_l 
vanishes at 0{h) if one sets oc = {—to, 0, 0, k). 

This leaves us with Uz, which we will determine by solving the remaining u = z,0 components 
of the stress tensor conservation law, to linear order in the metric fluctuation, and to leading order 
in gradients. First we compute 



TTZZ ^ _4 

-'-'■1st — s'l 



Oz) 



c 



Oz) 



(33) 



where h±_i_ = h^x + hyy. Then V^T^^ = yields 
while from V^T^° = we get 

K ^ - \K - (35) 

Specialize now to conformal theories (^ = 0, P = |e). If we stay away from the sound pole 
^ |A;^, then to leading order in derivatives we can ignore liH^ and the terms in (|34] ). and 



combine dM]) and dSSj) to get 



and 



Uz^^ATQ,i\Ko + \K.-\Kz). (36) 



' = ^ + I g2_ 1^2 (^0. - 1^00 - l^n) + 0{h\uh, kh). (37) 

t 3 ^ 

In what follows we will derive Kubo relations for the second-order hydrodynamic coefficients 
from the response of the stress tensor component T^^. From ([T) and Q, 

T^y = e [|f/(i)-f/(i)^ + li-Ky + K^v'^^'Ky)] + l{5^'h)i-Ky) + wy + 0{h'), (38) 

where the superscript (1) denotes a contribution of order h and e = e + 5^^^e + ■ ■ ■ . 



*Our notation 0{ujU±,kU±) assumes that the leading order in U± is not zero. 
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2.2 Ai 

In the hydrodynamic expansion, Ai is the shear tensor-squared coefficient: 



(39) 



where we will use the notation (• ■ • )ai to indicate that we are only showing terms of ■ • • that 
depend on Ai. We will also need 



±2) 



(40) 



and 



'XX I ^yy — 3(2^^-L ^zz) 3(^2 ^02^ 



xy — K'xy 



(41) 
(42) 



The Ai dependence (T^^)ai of (|38l) comes solely from the term at the order shown: 



^xyy^xx ~l~ ^yy) ~^ ^xz^yz 



Ai 



2(f/; - /^oJ) + (t/:. - + /i.,)(f/; - h' + V) 



(43) 



We now want to devise a simple Kubo-like formula for extracting Ai from some limit of a retarded 
3-point correlator. To that end, note that in the limit of x-independent sources, h^v = h^y{t), 
equation (l43l) simplifies to 



In momentum space this becomes 



Ai 



lKy{\h±^ 



hzz) ~l~ hxzh 



xz"'yz 



(44) 



-Ai / d'^qi / d'^q2 



X 



hxy{qi){\h±^{q2) - hzz{q2)) + Kz{qi)hyz{q2) 



(45) 



where the 4-momentum q = {to, 0, 0, k) is associated with the measurement point, and the 4- 
momenta qi = (wi, 0, 0, A;i), ^2 = (w2,0,0,A;2) are associated with each one of the h factors 
respectively. By (fTTI) . this produces contributions to G"*^'^^'^-'-, G"*^'^^'^^ and G^^^^^^^^ in the limit 
where all three spatial momenta vanish. For the purposes of our later AdS/CFT calculation, we 
will find that the least technically challenging case to calculate will be G^^'^'^l^^. So here we will 
focus on finding a Kubo-like formula for Ai in terms of G^^'^'^l^^. More Kubo-like formulas can be 
found in Appendix |Bl Consolidating identical terms in (fTTI) gives 

{T''{q))k = I d% I d%-^^6\q -q,- q^) -q,, -q2)Kz{qi)hyz{q2) 

+ other G's. (46) 
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In what follows we will suppress writing the momentum dependence of the 3-point correlator, with 
the understanding that the momenta —qi, —q2 are associated with the second pair and third pair of 
indices of the 3-point correlator respectively. We are now ready to compare (|46] ) to (|45] ) to extract 
Ai: 

lim d^,d^, lim G^J^I^^I?^- = -Ai + (Ai-independent terms). (47) 

oji— s-O fci— >0 

U12—I-0 fc2-i-0 

We leave the derivation of the A i -independent terms to Appendix |Al The result is: 

hm d^,d^, lim G^^l^^lf- = -Ai + r/rn. (48) 

2.3 A2 

To evaluate the A2 contribution to the stress tensor T^^, we need to compute 

Assuming, as we did before, that the fluctuations depend only on t, z coordinates, vorticity is 
given by 

/ -f/^,^ 

m ~ n'^ ~ |(9,f/, - djU,) = U -f/; I . (50) 



Substituting dSO]) and (1401) into (1491) we find 

(T^y),, ^ iA2 ( - U'^iU'^ - + hy,) - U'yiU'^ - + h..)^ . (51) 

Again, we are interested in finding a simple Kubo relation which is amenable to a straightforward 
AdS/CFT computation. We settle on using G^y\y^\^^ and we consider sources such that hy^ = 
hyz{t) and ho^ = ho^iz). After a substitution of U_l from (13T1) into (ISTT ). and a comparison with the 
relevant terms from (ITTl ) 

(T^^g)). = 1 1 d%-^^6\q -q,- q,) G^y^y^^'^ {q; -q^, -92) V(qi)/^o.(g2) 

+ other G's (52) 



we are led to 

lim dk^d^, lim = -iAs + ^tu, (53) 

where the A2-independent term ^t]Tyi is derived in Appendix lA] Other Kubo-type relations can be 
found in Appendix IE 
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2.4 As 

A3 is the coefficient of the square of the vorticity tensor in the hydrodynamic expansion of the 
stress tensor. To get A3, we rely on computing the same component T^^: 

(T-^),3 ^ (n-^),3. ^ iA3t/^f/; (54) 

Substituting (|3TI) we get 

{I ^jAs - -4A3fcifi;2 ^- , r, ,2\ r~ ; n ?2N • (5^) 

A simple way to extract A3 from terms up to second order in momenta in the 3-point correlator (at 
higher order in momenta we would need to expand the stress tensor to terms that include higher 
order hydro coefficients) is to take the static limit: coi = 102 = 0. Then, a comparison of (|55l) with 
the stress tensor response in terms of Green's functions ([TTI) . namely 

(T^^(q)), = I d% I rf^g2^54(g-qi-g2)G^'^'°^'l°nq;-gi,-g2)/^o.(qi)V(q2)+otherG's, 

(56) 

yields 

limdkA, lim = -iA3. (57) 

k2—>-0 uJ2^0 



We will use AdS/CFT to evaluate A3 from (|57l) . However, for the reader's convenience we give 
other Kubo formulae in Appendix |Bl 

3 Stress tensor correlators via AdS/CFT 

The generating functional for the stress tensor correlators of A/'=4 super Yang-Mills theory in the 
limit of large number of colors (N^ 3> 1) and at strong coupling is the on-shell five-dimensional 
gravitational action composed of the Einstein-Hilbert action, a cosmological constant term, the 
Gibbons-Hawking term, and holographic renormalization counterterms [.201 1271 122ll23l : 
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JdM JdM a — Z Jqj^ 

(58) 



^-AdS IJM 

where the values of the cosmological constant and of the volume counterterm parameter a are 

A = ^-^ — and a = — ^ '-, with d = A. (59) 

o r?2 D ' ^ ' 

The trace of the extrinsic curvature tensor can be expressed in terms of the induced metric on the 
boundary {gbdy)MN = guN ~ riMnN and the unit normal to the boundary n'^^ as 

K = {gMN - nMnN)V^'n^. (60) 
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In general, Ai asymptotes to an AdS space of radius -RAds- Since we are studying finite-temperature 
super Yang-Mills theory, the background Ai is five-dimensional AdS-Schwarzschild space [1241 : 

ds^ = + X-) + f{u) = l-u\ (61) 

Taking advantage of the fact that the AdS radius drops out of final results, it is convenient to set 
-^Ads = 2 and work in units 2tiT = 1 so that 

ds^ = + -YJT^- (62) 

Next we will evaluate the on-shell action, by expanding the metric in fluctuations Sqm n (where 
M, N = 0, 1, 2, 3, 5 = t, X, y, z, u) around the AdS-Schwarzschild background. In imaginary time, 
AdS-Schwarzschild is smooth and non-singular, and 3-point correlators are given by Witten dia- 
grams as depicted in Fig. la. In real time, as discussed in [|25l [T5l , the corresponding diagrams 
for retarded correlators G^aa live in the right-quadrant of the Penrose diagram due to causality. 
The lines in this diagram represent advanced and retarded bulk-to-boundary propagators. These 
propagators are identified as solutions 5(7^ to the linearized equations of motion, which approach 
prescribed values at the boundary n = and which reduce (up to gauge terms) to purely incom- 
ing/outgoing waves: e~*'^*(l — n)^*"^/^. 



3.1 Two-point stress tensor correlators and comparison with hydrodynamic 
expansion 

First we will warm up with 2-point retarded correlators. We recall that the gravity fluctuations 
are taken to be independent of x, y coordinates, and are slowly varying functions of t, z. Given 
the symmetry of the problem, the fluctuations can be classified according to their transformations 
under rotations about the z axis as: 5*0(2) tensors {Sg^y and Sg^x — ^9yy)^ vectors {dg^x-, Sgzx and 
Sgoy, Sgzy), and scalars (all others). 

For completeness we list in Appendix O the bulk metric fluctuations in momentum space as a 
series expansion mco,k, and in terms of the boundary metric fluctuations, h"^. The retarded bulk- 
to-boundary graviton propagators are easily obtained from the expressions given in Appendix O 
by differentiating with respect to the boundary fields. We work in the gauge 

SgM5 = 0. (63) 

The tensor modes 5g!l. and 5g^ — 5g^ have propagators which do not have singularities for 
u,k <ti 1 lfT9l . The vector (or shear) modes, e.g. 5gQ, 5g^, have poles typical of diffusion, and the 
scalar mode propagators have a sound pole [|26l . We will refer to the vector metric fluctuations as 
shear modes, and to the scalar metric fluctuations as sound modes. 

The computation of the stress tensor 1 -point function is reviewed in Appendix |El 10: 

o/loo 2'7r^ 



^Restoring the units, the energy density and pressure equal e = g" — , P 
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gives the energy density of the finite temperature strongly coupled super Yang-Mills theory and 



{T% = 2^ = S'^^ ^ tj = 1,2,3 = x,y,z, (65) 

gives its pressure. 

The two-point retarded correlators are computed using the quadratic gravitational vertex re- 
viewed in Appendix E.2. Since we are interested in ra correlators, one of the bulk-to-boundary 
propagators is retarded (the one associated with the largest time point on the boundary), and the 
other propagator is advanced. However, we should remember that the bulk-to-boundary propaga- 
tors in momentum space are defined as retarded or advanced relative to the momentum conjugate to 
the boundary space-time point. 4-momentum conservation gives q = —qi, where q is the momen- 
tum associated with the point which has the largest time. This effectively transforms the advanced 
bulk-to-boundary propagator Ga.riqi, u) into a retarded propagator Gva.iq, u). 

From the decoupled fluctuation 5gy one recoversl^ the stress tensor retarded 2-point function 

m 



_fS_ 

+ ^ + ■■■ (66) 



27^2 2V 
The result derived in the hydrodynamic limit from solving dS]) and dH) is 

Glto = ¥ - 'Voo + vruco' - + k') + ■ ■ ■ (67) 

where dots denote terms which are higher order in gradients (and hydrodynamic expansion coeffi- 
cients). By identifying (|66l) and (|67l) one obtains a handful of hydro coefficients, of first and second 
ordeiO: 

iV2(2-ln2) 

^ = ^' ^ = ^' ^^n= ' g 2 ■ (68) 

From (|64] ) and (|68T ). the diffusion constant of the shear modes in the strongly coupled J\f=A super 
Yang-Mills plasma is 

ri 1 

D„ = — '-^ = -. (69) 

^ e + P 2 ^ 

This is the maximal set of hydrodynamic coefficients which can be determined from 2-point stress 
tensor correlators in the AdS-Schwarzschild background [2]. On the other hand, having found 
e, P, 1], tyi and k, we should be able to match all the other 2-point correlators obtained via AdS/CFT 
with their hydrodynamic expressions. 

To see how this plays out we consider the 2-point functions of shear modes. We perform an 
expansion assuming a shear dispersion relation w ~ /c^ ~ ^ 1, where A is a small expansion 



^For the sake of brevity, we suppress energy-momentum conservation delta functions. 
Restoring the units, this reads r/ = — | — , k = — , r/rn = — - — . 
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parameter, and extend the results of Policastro, Son and Starinets [fT9l to second order in A: 



G 



Ox I Ox 
AdS 



6hn,.6h 



2%^ 



k"^ + 6iuj 



Ox\xz 
AdS 



6^S 



Ar2 

c 

2%2 



22(-iu; + ip; 



A;2(2a;2(ln2-l)-2ia;fc2 + fc^) 



+ ■ 



G 



xz\xz 
AdS 



ShxzShxz 



2%2 



1 

23 



+ ■ 



(70) 



As a result of this expansion, the shear poles are visible. The presence of the higher order poles in 
the 2-point correlators is an artifact of the expansion and indicates that the location oj = —i^k"^ of 
the shear pole is shifted by higher order terms 
This can be contrasted with 



G' 



Ox I Ox 
hydro 



y^Ox\xz 
^ hydro 



^X2|X2 

hydro 



—ioj + Dj^k'^ 
rjuk 



+ 



ujk'^{-{r]Tu - + jiKD^k"^) 



Dr,k^y 



+ 



-iuj + Drjk"^ 

,2 



+ 



rjtu 



-lUJ 



+ 



(71) 



There is an apparent mismatch in the hydrodynamic expansion and the AdS calculation in Q^^\'^^ ^ 
namely the coefficients wA;'^ and k^ which multiply the second order shear pole disagree. This 
puzzle is resolved by noting that the expansion in small w, k assuming a shear dispersion relation 
has the rather unwanted effect of mixing hydrodynamic coefficients of different order in each term 
of the series expansion: e.g. e and 77 at zeroth order in A; r^rn — |k and D^^k at second order in A; 
etc. So a contamination with third-order hydro coeffcients is to be expected, as forecast by the term 
proportional with which contains the product of two second order hydro coefficients. This 
is further elucidated by a comparison with what happens when instead expanding the correlators 
assuming a; ~ /c ~ A <^ 1. Then, the AdS result is 



a 



a 



G 



Ox|Ox 


Ar2 

c 


'3 


AdS ~ 


2%2 


23 


Ox|x2 


Ar2 

c 




AdS ~ 


2%2 




XZ X2 


Ar2 

c 


' 1 


AdS ~ 


2%2 


23 



+ 



22a; 
k k(2u?{\ 



22 



+ 



23a;2 
In 2) - P) 



+ ■ 



'2^ 



+ ■ ■ 

(2a;2(i_i^2)-P)7r2Ar2T^ 
23 



+ 



(72) 



'^The shear pole, as extracted from second order hydro, is [)2j: w : 
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while the hydro expansion gives 

e - I \ ^ \ 

UJ 

—irjk H h . . . 

u 

|e — irjLO + r]TiiCO^ — — rjDjjk'^ + ■ ■ ■ (73) 

The UJ poles are a remnant of expanding the shear pole, but now the two sets of 2-point correlators 
are perfectly matched, and at each term in the expansion there is no mixing between different 
order hydro coefficients. One may extract e, r] and the combination r^r^ — with the same result 
as before. In this respect, the tensor mode 2-point correlator l^s^ enables the identification of a 
larger set, since r]Tn and k can be obtained separately. 

In appendix |D] we list the 2-point functions of sound modes, computed using a sound mode 
dispersion expansion a;~A;~A^lto second order in A, extending the results of Policastro, 
Son and Starinets [[26l . The higher-order poles which appear in appendix iDl are also an artifact 
of expanding the higher-order attenuation terms which are present in the sound pole. There is no 
mixing of the different hydro coefficients at each term in the expansion, and the hydro and AdS 
results match, provided that one takes into account that the AdS computation yields tensor density 
correlators (2{6'^S) / {5h^,y5hpa-) = 5{y/^Tpy) / 5hpa) and not tensor correlators (5{T^'') /5hpa-). 
The difference is a contact term, namely 

= ^rds"' + V'^'T^", (74) 
which arises from differentiating the volume factor with respect to the metric fluctuations. 



^Ox|Ox 
hydro 

^Ox\xz 
hydro 



hydro 



3.2 Ai , A2 and A3 via AdS/CFT 

Due to the complexity of evaluating generic stress tensor 3-point functions, in this section we 
contend ourselves with computing lim/ci-^^o G^^g^'^^, limtJi->.o and lim.^^o G^^^^'^^ We 

leave a more complete computation of 3-point stress tensor correlators at arbitrary (albeit small) 
values of 4-momenta for future work. 

Note that we are using 3-point correlators where one pair of legs, xy (which is associated 
with the largest-time boundary point) corresponds to tensor fluctuations, and the other leg pairs 
correspond to 0(2) vector (shear mode) fluctuations: xz,xO and yz,yO. Because of the mixing 
between Sg±o and 5g±z, each of the corresponding bulk-to-boundary propagators is a 2 x 2 matrix. 
So, for example, the Witten diagram for G^^/ls^^l^'^ has four terms, depicted schematically in Fig. 2. 

Moreover, each 2x2 shear-mode propagator contains a pure gauge (diffeomorphism) contri- 
bution whose behavior at the horizon is not in the form of an incoming/outgoing wave. This gauge 
artifact can (and does) cause technical difficulties in integrating the location of the vertex in the 
Witten diagram up to the horizon^. We will avoid this issue by computing correlators for which 
the diffeomorphism terms drop out, that is, we will be focus on calculations where we can easily 
work with gauge-invariant modes. 

^The problematic behavior is related to the fact that the diffeomorphism parameters (as constrained by the gauge 
condition 6gM5 = 0) are non-analytic functions at the horizon (see equation (lllll) ). 
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(a) (b) 




Figure 2: Witten diagrams for the retarded 3-point correlator G^^^l^^^l^^ with the boundary point x 
having the largest time; Xi and X2 can have any time order. 



As an example, the low momentum solution to the equation of motion for and 6gl which 
gives the retarded bulk-to-boundary propagators is0 

+ (l-^l)\n{l + u)) + ■■■]-£),- (75) 

Z J bJ 

where 

= ^■■■yK^-^K) 

~ _ ( e k\2A2 + e) \ 

— i^-yi (76) 

The Di terms in (|76l ) represent a gauge mode associated with shifts in the x direction. One possible 
way to trivially isolate the gauge mode is to take the A; — )■ limit. Then the the shear modes 



"'Here, since we are interested in sending the spatial momentum to first, we give the shear mode propagators 
assuming a sound mode dispersion relation. In Appendix ICl we assumed a diffusion dispersion relation to highlight 
the presence of the diffusion pole. 
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decouple. 5(?g(A;=0) is equal to the boundary value up to 0{uj^) corrections, and 5g^{k=0) is a 
purely incoming wave at the horizon to the same 0{uj^) order: 



6g:ik=0) = h:ik=0) X ( 1 - ^'(^ ^""ll "'^ + • • • ) (1 - n)-/^ (l - ^^ ln(l + u) 



+ a;2(_lLi2(i^) + lln2(l+«) + (l-l^^)ln(l + «)) + ---^. (77) 

As a result, only Fig. 2a contributes to G'Xds^'^^l^i"'-'' ^2=0) [a consequence of 0(3) rotation 
invariance when A; = 0], and we expect no subtleties coming from the near horizon region. 

Except for the specific form of the cubic gravitational vertex, the evaluation of the 3-point stress 
tensor correlator is no different than evaluating a causal scalar 3-point correlator with a derivative 
cubic vertex in the AdS-Schwarzschild bulk0. Consider this scalar example, with cubic interaction 
Xvf V~99^^^'^M<pdN(l), for the sake of simplicity of discussion, and examine the small frequency 
(with k's zero) behavior of this contribution to the action S. As we'll see, this is dominated by the 
near-horizon (u 1) contribution 

v/=^^7">a„(/.a„</.. (78) 

The small frequency behavior of the retarded scalar bulk-to-boundary propagator is (1 — m)^"^/^. 
Therefore, computing the small frequency behavior of the integrand in (|78T ) gives a contribution to 
the causal retarded correlator which in the near-horizon region is proportional to 

du{l - u)-'^'^'+^''^-\uiU2 - (Wl + ^2)2) (79) 



where the terms quadratic in w's arise from taking the u derivatives. In arriving at (1791 ) we took the 
momenta flowing through the a legs of the raa correlator to be coi and uj2. By energy conservation, 
the momentum flowing through the r leg is cu = —{ui + U2). The potential divergence at the 
horizon is regularized by the fact that the causal 3-point function of interest corresponds to cui — )■ 
cui — ie, UJ2 ^ UJ2 — and a; = — (cui + 102) — (cui + + 2^elll Integrating (|79l) from u = 
to 1 yields 

{UJ1UJ2 - {uji + ws)^) (80) 

wi + 002 

to leading order in w's. Counting powers of frequency as A ~ ~ a;2, the near-horizon divergence 
of the integrand has enhanced the naively O(A^) order of the integrand to an 0(A) contribution to 
the 3-point function. 

A similar story emerges for the G^jg^'^^(A;i = 0, A;2 = 0) correlator. There, the integrand, as 
given by the expansion of the cubic gravitational action, behaves near the horizon as 

2^ V (1-u) (1-u) J 



''if the scalar supergravity fields have a derivative coupling of the type \J—gg^^^ dM4>dN4>, then, for on-shell 
fields this reduces to (to^/2) \/—g4>^ plus a total derivative term. However, in the text, we refrain from performing 



IM 

the integration by parts and merely comment on the behavior of the integrand near the horizon 
See the discussion in sec II.B of [25 J . and see I27n28il . 
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where at least one of the tu factors arise from derivatives on the wave factor (1 — u) «('^i+'^2){^ 
Regularizing the wave factors as before, and performing the integral over the black hole bulk up to 
the horizon, gives a finite contribution to the G'^jg^'^^(/i;i=0, /c2=0) correlator: 



lim G 

fc2-i>0 



xy\yz\xz 
AdS 



Ar2 

c 

2%2 



1 

23 



.U1+U2 (a;ia;2 + + i^2)(ln2 - 1) 



fcl-5-0 



fc2-5-0 

to get Ai: 



Ai 



(82) 



The last step is to identify the AdS/CFT correlator with the corresponding hydro expansions 

lim Ghycto '""^ = ¥~ ^Vi^i + ^2) + V^niuJi + ^2 + ^1^2) - hf^i^f + (^2) - >^i^i^2 + ■■■ (83) 



(84) 



For the remaining A2, A3 coefficients we will need to take the limit w — )■ of the shear 
modes. From Appendix O we get that Sg'-^{uj=0) is equal to its boundary value to order 0{k^), 
and (5(^0 (a;=0) decouples: 



-k'^u{l-u) 



(85) 



Thus, once again, the graviton shear modes decouple, and the correlators needed to evaluate A2 and 
A3 are given by a single Witten-type diagram, with the causal graviton propagators being diagonal 
in both the tensor and shear mode sectors. The integral over the cubic action is straightforward and 
gives: 



lim G 

a;i->0 
0)2— !>0 



xy\0y\0x 
AdS 



Comparison with the hydro expansion 



lim G 

W2— >-o 



xy\0y\0x 
hydro 



c 



\ + {kl + kl) + 



\t + \K{kl + kl) - \\^kik2 + 



(86) 



(87) 



yields 



A3 



0. 



(88) 



'^As a result, finding the integrand in (ISTt only requires knowing the solution Sg'^'^^ to the linearized equations of 
motion to 0(a;J, a;|), as given in ( |79] |. 

''^ Given the discussion preceding (|74| |. readers may wonder why we are identifying tensor density correlators 
(GAds) with tensor correlators (Ghydro)- There is no difference for the particular components we are interested. In 
general, 

'(^, ^) - G';:}f^\x, y, z) + Gl^^fix, y)g<5{x - z) + G^Js ^(^, ^)9''^Six - y) 
+GZsi^)i9"'9^^ + + 9''^9nS{x - y)d{x - z). 
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Similarly, we compute 



which can be identified with 



limG** = ^..fc. + ... (89) 



lim = (-iA, + lr,rn)i^ife + • • ■ (90) 



fci^O 



to 



yielcS 



Ar2 

^2 = -7^ In 2. (91) 

This concludes our derivation of the second order hydrodynamic coefficients Ai, A2 and A3 from 
3-point stress tensor correlators^. 
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r]Tu contribution to Ghydro ^hydro 

In this section we are interested in computing the r/rn contribution to T^y to second order in gra- 
dients, and in the presence of a gravitational perturbation which includes only the shear modes: 
hxo, hxz, hyo and hyz. We begin with 



(n^^),^ = r]Tu (^^UVa^'-"^ + I V ■ f/ a^'^^ (92) 



'^See f^^l for arguments regarding a universal relationship between Ai, A2 and rjTYi- 

'^In addition to the computation detailed in the main text, we have performed a separate check of (IT4t using a 
different set of gauge-invariant modes e.g. choosing the boundary fields such that h% — —h^ or /ig = —^k'^{h^ + 
hv) + ■■■ ^ {-1 + P){h% + /ip + • • • = -^{ik + kuj\n2){h% + hy) + ■ ■ ■ where ellipsis denotes that 
we have required that the gauge component of the bulk-to-boundary propagators vanishes to order 0{X^) in a small 
uj ^ k ^ X 1 expansion. For example, solving the hydrodynamic equations we compute 

G™'°^l"^(a.i-a.2=0) = -le + nikj + fc|) - ^X^k.k, + ■■■, 

Gh;d?o'"''"'"'"'(^i=fc2=0) = |e-4^ryK+c.2)+4ryr^(a.?+^| + ^la.2)-2«2(a.?+c.|)-4Al..lC.2 + ••• , 

G^;,7/I^-I°^(fci=0, c.2=0) = i-lX, + ^ru)io,k, + ■■■. 
The corresponding AdS 3-point functions are 

'-'AdS lWi-a;2-Uj — -26;;^ + js^j^V'^i + fcjj H , 

xx+yy-2zz\xy\xy,, , _„^ _ N'^ C, , 4- , , ^ Af,^(ln2-1) ,2 , , ,2 , , , , , N , 

and 

Gxx~yy\zx\^x / -, ^ ^\ A^^ , 

AdS (fci=0,W2=0) 2j^a;i/c2 + ••• . 
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where 

<[/Vcr^^> = A'^^A^'^t/^V^fTp^ - iA'^^A^'^t/^V^cTp,. (93) 

Given the set of metric perturbations considered, cr°^' = 0{h'^). The non-vanishing components of 
cr*-^ , to linear order in the shear metric fluctuations, are a^^ and a'^^ given in (|40l) (where we recall 
that the shear tensor a is symmetric). So for (n^^)^T^n o'^^Y '^he first term in (|92] ) is non- vanishing. 
We will also need the projectors A^*^ to order 0{h'^): 



/o 




Uy 


\ 




1 







Uy 





1 


-h 






~hyz 


1 / 



+ 0(/i2) (94) 



where U^, Uy are given in (|3T| ). Given that A"^^ vanishes to 0{h^), only the first term in (|93T ) for 
IJ, = x,iy = y contributes to terms quadratic in the metric fluctuations: 

do^T^J - m^^yz + ny^.z) + (A^^^oa,, + A^^^ofT,.). (95) 
The last piece we must compute is a^y to second order in the fluctuations 



g = /ixMO-(i) + h,ya^'^ + (96) 



where the truncation of (H) to second order in fluctuations. Using that the second order 

expansion of the Christoffel symbols is 



r(2)A^ = ~g'^h,„T^^)l^, (97) 



(2') 

the shear tensor axy evaluates to 



a. 



(2) 
xy 



do{UxUy)+oihe,rh'^. (98) 



Substituting everything into (|92l ) we arrive at 

(n-^),,, = r/rn [dHU^Uy) - \{Kx - K^W^ - h',^ + Ky) - \{Ky - h',y){U', - + hx) 



hxziU'y - Ky + Ky) - hy,{U'^ - h'o^ + Kx) + other h's^ . 



(99) 



Recalling that Ui_ velocities vanish for z-independent metric perturbations (that is, Ui_ vanishes 
at zero spatial momentum), we can now extract the desired r^rn contribution to the correlators 

G';t'^^(A:i=0, ^2=0) and G^''{k,=Q, u,=Q): 

G'hydro '^""(^1=0' ^2=0) = riTii{uJiUJ2 + ul + ul) ^ h r^Tn "independent terms (100) 

^hydro '°''(^i=0' 1^2=0) = \riTj1k2UJ1 H h r^rn-independent terms. (101) 

More general explicit computations show that there are no other second-order hydro coefficients 
contributions to (|48l |53l |57] ) besides the ones we explicitly went over in Section 2 and in this 
appendix. 
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B More Kubo formulae 



Using the method outlined in Section 2 one can derive a variety of Kubo relations. We give 
below a few more to supplement (|48T ). (|53T ) and (l57l) . To simplify notation we suppress writing 
the momentum dependence of a retarded raa 3-point correlator G^^^^^^^^{q = qi + q2] — Qi, — ^2), 
with the understanding that the momenta —qi,—q2 are associated with the second pair and third 
pair of indices of the 3-point correlator respectively. 



B.l Ai 



aJ2->'0 fc2-s>0 

lim 9^,9^2 lim = -|Ai - |k + r/rn 

U}2^0 fc2-S>0 



lim 9^1(9,^2 lim Ghvdro = -Ai + ?7rn. (102) 



(^2^0 fc2-s>0 

Noting that 0(3) rotational symmetry is restored in the limit of vanishing spatial momenta, we can 
use (11021) together with 0(3) covarianceto compute other correlators. E.g. G^^^^^^^y{ki=0, ^2=0) = 
Qxx\xz\xz(^)^^^Q^ A;2=0) = G^^I^'^l^'^(A;i=0, A;2=0) and G^^I^'^I^J'(A;i=0, A;2=0) = G™l^^l^^(A;i=0, ^2=0). 
Rotational symmetry can also be used to check the last Kubo relation in (|102|) : a 7r/4 rotation in 
the {y,z) plane can be used to show that G2^^l^?/l^^(A;i=0, A;2=0) = i[G^fl^^l^^(fci=0, A;2=0) - 
Gyy\^y\^y(^ki=0,k2=0)]. 



B.2 A2 

The order of limits (t^i — )■ and k2 — ?■ 0) is important here (due to the presence of a sound pole in 
the correlators). On the one hand, 



lim dk^ lim d^^ Urn O^ = -A2 - + 2r]Tn 

k2^0 LUi—>-0 fci— >-0 ■' 
012—^0 



Hm dk^ lim lim 0;^^'''°'' = -^Aa - n + r]Tu (103) 



a;2— ^0 
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but a different order of limits yields 



lim duj^ lim dk^ lim G^d^ '""^ = -|A2 - k^. + h^u + ^ 
wi-^o k2^o fci^o ^ ^ 2 3' 2e 

0^2— >-0 

lim lim lim Ghvdro = 1-^2 + ki^ - Iv^u - — 

uJi^O k2^0 fci^O 2 3' g 

aJ2->-0 

lim a., lim lim = ^A^ + - ir^m - ^. (104) 

a;2->-0 



B.3 A3 

nma.,a., nmG;tr^ = |A3, 

lim 9,^5,, lim = iA3. (105) 

A:2-!>0 ti;2-!>0 

In the relation (l53l) given in the main text, one gets the same answer independent of the order of 
limits wi — )■ and k2 ^ 0. 



C Retarded bulk-to-boundary propagators 

The metric fluctuations, solutions to the linearized Einstein equations, are in general coupled. We 
work in the gauge 

6gM5 = 0, M = 0,l,2,3,5 = t,x,y,z,u. (106) 

We define 

59^' ^ -g'^'Sgp^ (107) 

The fluctuations considered in this paper are independent of x,y coordinates. Relative to the 
rotation group 50(2) about the z axis, one distinguishes tensor fluctuations {5gxy] 5gxx — Sgyy), 
vector fluctuations (5 got, Sg^z ; and 5gyo, 5gyz, and scalar fluctuations (5 goo, 5gzz, Sgoz, and Sgxx + 
5gyy). The tensor fluctuations are completely decoupled, while the vector and scalar fluctuations 
all mix within their respective sectors for generic uj,k. For the vector, 5gxo mixes with 5gxz and 
6gyo mixes with 6gyz. 

Since we are interested in the hydrodynamic regime, where the fields vary slowly with t, z, 
the equations of motion are solved in the bulk perturbatively in tu, k. The bulk-to-boundary causal 
propagator of coupled functuations will contain one term which behaves like an incoming/outgoing 
wave at the horizon, in addition to terms which are diffeomorphism terms. The existence of these 
diffeomorphism terms is inferred by solving the equations of motion near the horizon to leading 
order, substituting an ansatz of the type F{u){l — uY as in [|26il . The values of r = ±iuj/2 
correspond to the incoming/outgoing waves; the other possible values of the exponent (e.g. r = 
0,-1/2) correspond to the diffeomorphism terms. Lastly, we require that the bulk fields approach 

prescribed values at the m = boundary: 6g^j h". 
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Wefin43 

= C,{1- ur-/' (^1 - ^1 ln(l + u) + co'i—u/-^) + 1 ln^(l + u) 



Sg 



Sgo 
Sgo 
5gl 



+ ln(l + u){l-\ ln2)) + P(i ln(l + u) + - u)) + . . . 
Z o o 

+ C,{1- u)-^-l^ {\ - .| ln(l + + c.2(-iLi2(^^) + ^ ln^(l + u) 

+ ln(l+M)(l-^ln2))-A;2ln(l + u) + ...^ - D^^/f 

= - (l - z| ln(l + n) + a;2(-iLi2(i^) + ^ ln^(l + n) 

+ ln(l + - i ln2)) + ln(l + m) + ^(1 -«)) + .. 
Z o o 

- - u)--/' (l - z| ln(l + u) + co\-^-U,{^-^) + 1 111^(1 + u) 
+ ln(l + M)(l-iln2))-A;2ln(l + M) + ...^ - D^^/f 

5gl = C,il-u)-'^/^(^-2 + tu\n{l + u)+u;\U2{^^)-^\n''il + u) 

- (2 - In 2) ln(l + m)) + k\-'^ ln(l + m) + ^(1 - m)) + . . . 



2k I— 
D3— + £'5(- V/ + 2A;^ arcsin(M)) 

CO 



Clil - u)-''^'^{-ikf + uk{l -u){u+ ^(1 + n)(ln(l + u) - 21n2) + . . . ^ 



+ 1^3 + -DoT^ D^ujk{u\/f + arcsin(u)) 



C3(l - - ,/ + |/ In + (u; + ^A;2)«(l _«) + ... ^ + 

CH - ur^/^k - ^/ + |/ In + + tk^)u{l - tx) + . . . ^ + 

C3(l - u)-'^'^ (^-2 + iu\n + A;' + . . . ^ - Z^i 



5c/^ = C4(l - - 2 + In + ^ + . . ^ - D 



'^We recall that, for simplicity of notation, we work in units where 2'kT = 1. Alternatively, one should think of lo 
and k as energy and momentum made dimensionless by division with 2'kT. 
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+ l\n^{l + u) + {l-^-^)\n{l+u))-kHn{l + u) + ...] (108) 
o 2 y 

where the coefficients Ci, . . . C^, Dq, . . . , D^, are given in terms of the boundary fields as fol- 
low ^^l : 

/ . 2uk^ k^A \ 

' ~ \3uj^-k^ \3a;2_p)2 + 24(3a;2_p)3 +■■ J^o 

^ V 3w2_p +'(3^2_p)2 24(3u;2 _ p)3 + • • J - 

/ 2ujk . 4a;2^3 \ 
^ V 3a;2_^2 +^(3^2 _ ^2)2 " I2(3a;2 _ p)3 + • • • J ^0 

A = 9uj%7r^ - 6 In^ 2) - 6a;^A;2(7r2 + 24 In 2 - 6 In^ 2) + ^^/^^(Tr^ + 48 In 2 - 6 In^ 2) - 32/^*^ 

_ / i/t a;A;((zw + ip)ln2-A;2) \ u; 

= UHc^ + I^^) 4Ha; + iP)2 + 

_ / .fe c.M(zc. + ifc2)ln2-A:2) \ . , , 

^ " ^2Hu; + iA:2) 4Ha; + ifc2)2 +--Jl^o+^^J 

^ (l. -^^-^;;^-^^^ ->- (109) 



/ 3a;2 2a;fc^ 2a;2^4^3^2^2 - ln2) + fc2(in2 - 4)) \ ^ 

° ~ V3a;2 _ p \3u;2 _p)2 + (3^2 _ p)3 

u;2 .ujk\u^-k^) w2p(^2_p)^3^2^2-ln2) + A;2(ln2-4)) 



_3u;2_p (3^2_p)2 (3w2_p)3 



UJ 



2 O, ,31,2 



2w3A;2 2w^A;2(3u;2(2 - In 2) + A;2(ln 2-4)) 



^ ' ~ 3a;2_p ^\3w2_p)2 (3a;2 _ p)3 +•••' - 

2wfc . 4^2^^ 4a;3A;3(3a;2(2-ln2) + fc2(ln2-4)) 
^ ' ~ 3tj2_p +^(3^2_p)2 (3a;2 _ p)3 +••• ' " 

_ A; P(u;ln2 + (V2)A;2) \ , 

^1 - 2Hu; + iP)2 +--J^o 

fc2 A;2a;(a;ln2 + (z/2)fc2) 

' 2(-«a; + iA;2) 2(-za; + iF)2 



'^The appearance of higher order sound and diffusion poles is only an artefact of expanding in small w, fc of simple 
poles which have additional attenuation pieces. For example, C3 can be repackaged as fc'^ (1 — (i /2)u) In 2 + . . . )/ (fc^ — 
iuj{2-k^) + ...). 
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/ fc^ Pa;(a;ln2 + (z/2)fc^) \ 

_ / 3a; A; . ?>u^e ujk^{9uj^ In 2 - 3a;^fc^ In 2 + 2k^) \ . 

^ ~ V 2(3a;2_p) +*(3^2_p)2 + (3^2 _ ^2)3 o 

^ V2(3w2_p) +^ 2(3^2 _p)2 + 2(3^2 _ ^2)3 +••• 

f tuk _. 3u^k Lo^k{9ujHn2-3uj'^kHn2 + 2k'^) \ 

V2(3w2_^2) -'(3^2_p)2 (3^2 _ ^2)3 ^ 

/ 3a;2 , 6w3^2 2uj^k\9uHn2 - 3u^k'^ln2 + 2k'^) \ . 
^ V3a;2-fc2 ~\3w2_^2)2 (3a;2 _ p)3 o 

_ ( k"^ _. 2uk^ 2k\9uj\l-ln2) + 3uj^k\ln2-2)-k^) \ . 
^ ~ V3a;2_p ~\3a;2_^2)2 + 3(3a;2 _ p)3 o 

c^2 . a;fc2^^2 _ ^ fc2(a;2 _ fc2)^9^4(^ - ln2) + 3a;2p(in2 - 2) - A;^) 



,3a;2_^2 (3^2_p)2 3(3a;2 - A;2)3 



2 2a;3^2 2uj^e{9u\l - In 2) + 3^2^(1^ 2 - 2) - A;^ 



^ ' ~ 3u;2_p +*(3^2_p)2 3(3^2 _ p)3 +•••' 

2u;A; . 4^2^ _ 4wA;3(9w4^1 - ln2) + 3u;2A;2(ln2 - 2) - A;^) 
^ ' ~ 3u;2-A;2 ^ \3u^ - k^y 3(3a;2 _ p)3 +---)o 
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The diffeomorphism parameters which give rise to the diffeomorphism terms in the bulk-to-boundary 
propagators are: 

c f J , f -Dof , /7 -Di -^2 .D3 .D5A;arcsin(n) 



2uu ' ku' ku^ uu u ' Wa/J/ 

- ~ ii£'5A;arcsin(M),L)5V^M^ (111) 



That is, the "D" terms in (fTOSl) are given by (5^i^)diff = V^v^^^ + V^^^^v- The Do, ^1, ^2, ^^3 
terms in (II 111) generate boundary diffeomorphisms. The term generates an infinitesimal scale 
transformation in the vicinity of the n = boundary, plus additional boundary diffeomorphisms 
to preserve the gauge condition 5gMb = 0. Note that these diffeomorphism parameters induce 
singular gauge transformations at the horizon. This is a consequence of working in the gauge 
Sgub = 0. 
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D Sound mode 2-point correlators 



From the coupled fluctuations Sg^, Sg^, Sg^, Sg^, Sg^ (where for simplicity of notation we lumped 
together the sound modes with one more fluctuation, namely 6g^ — dg^ which can actually be 
decoupled) and from the form of the on-shell quadratic gravitational action (see Appendix E) we 
get 



-4- 



Nl_ 
2%2 



2 r 



— i- 



23(3a;2 - k^) 



(3a;2 - k^y 



+ 



1 



a;'(18a;*'(l -ln2) +a;V(241n2 - 39) + a;^A;*(28 - 12 In 2) 



2 7.4/ 



G 



xx\yy 
AdS 



2(36^2 _ ^2)3 

+ A;^(21n2-7)) ) + ... 



= -4 



^hxx^hyy 



2%2 



+ 



23(3a;2-A;2) 2(3a;2 - A;2)2 



9a;*^(ln2- 1) +a;V(15-31n2) -cc;4r(8 + 31n2) 



^00,00 _ . 

<^AdS - -4 

iV2 



2(3a;2 - k^f 
+ (3 + ln2)u;V -A;«) ) + 
5^S 



6^h 



00 



2%2 



3{5e - 3cc;2 

23(3a;2 - k^) "(3a;2-A;2)2 



— 



k\9uj\l - ln2) + 3a;2A;2(ln2 - 2) - A;^) 
^ (3a;2 - A;2)3 



G 



zz\zz 
AdS 



= -4 



6^S 



2%2 



3cc;^ 



(7a;^ - fc^) 

23(3a;2_^2) ^(3^2_^2)2 



cc;^(9cc;^(l - ln2) + 3cc;'fc'(ln2 - 2) - fc^)7r' 



(3C^2 _ ^2)c 



00|0^ 
AdS 



-2- 



5hnnSh, 



OO'Jil'Oz 



iV2 



2%2 



3(^A; 



_2(3cj2 _ A;2 



3a;2A:3 



(3a;2 - F)^ 



^ ^^A;3(9cu4(l - ln2) + 3a;2^2Qj^2 - 2) - k^) ^ 



(36^2 - A;2)3 
25 



ShooSh^z 
2%2 



— i- 



23(3a;2 - (3a;2 - A;2)2 



^k\9uo\l - ln2) + 3uj^k\\n2 - 2) - k^) 



(3C^2 _ ^2)3 



+ . . . 



5hozSh 



■Oz 



Nt 



(9^'2 ^ ^.2) ^ 3^,3^.: 



2%2 [23(3a;2 - A;2) (3cj2 _ ^2)2 
a;2A;2(9a;^(l - ln2) + ?,uj^k^i\Ti2 - 2) - A;^) 



(3a;2 - k^f 



+ ... 



-2- 



6hoz6h, 



Nt 



uk 



Su^k 



— i- 



2%2 [2(3^2 _p) "(3^2 _ ^2)2 

kuj^{'9uj\l - ln2) + 3(ln2 - 2)a;^A:^ - k^) 



(3a;2 - A;2)3 



+ 



2%2 



23(3a;2_^2) 2(3cj2_p)2 



- k^){9uj\l - ln2) + 3a;^A;^(ln2 - 2) - k^) 
2(3a;2 - A;2)3 



+ 



6^S 
) 

'6hnJh, 



QzOlOxx 



2%2 



+ ^ 



.3(a;2 - k^)kLJ^ 



2(3a;2-fc2) 2(3a;2 - fc2)s 



- e){9u\l - ln2) + 3(ln2 - 2)uj^k^ - k^ 



2(3a;2 - A;2)3 



+ 



-4 



{00^ + k^) 3{oo^ - e)oo^ 



2%2 L23(3a;2 _ ^2) 2(3cj2 - A;2)2 
u^jk^ - u;^){9u;\l - ln2) + 3(ln2 - 2)a;2A;2 - k^) 

2(3^^2 _ 



+ 



26 



These 2-point functions exhibit, as expected, the sound pole0 expanded in small uj,k: to 



E On-shell gravitational action 

For completeness, in this appendix we give the linear, quadratic and cubic gravitational vertex. 
This corresponds to expanding the gravitational action (|58T ) in terms linear, quadratic and cubic in 
the linearized on-shell fluctuations: 



s = s + 5(^)5 + 5(2)5 + 5(3)5 + . . . 



(113) 



E.l First order action 

To linear order in the fluctuations (|58T ) is a boundary term0: 



5(1) 



\ JM JdM / ^ JdM L \ 



(114) 



where jJi.y, p = 0, 1, 2, 3, and all indices are raised and lowered with the backgroud metric, e.g. 
Sg'^'^ = Sgp^g^'^g'^'^. As expected, the Gibbons-Hawking term has contributed to the cancellation 
of the terms linear in derivatives of the metric fluctuations at the boundary. Also, the leading order 
divergence (proportional to where the boundary value is of the radial coordinate ub — )■ 0) 

from the Einstein-Hilbert plus cosmological constant action and from the Gibbons-Hawking term 
is canceled by the boundary volume counterterm. 

In the black hole background, the on-shell action linear in the boundary fields is finite [|26]j^ : 



j=l,2,3 



(115) 



where we have introduced the shorthand notation 



u=0 



(116) 



u=0 



For example, the one-point function of the stress tensor is obtained by differentiating with 

,01 

Oln=0- 



respect to = 5g^ 



,55 



hY,=0 



3Nl 

277r2 



(117) 



'^The sound pole complete expression is 1 19] Vsk — 2{t+p) (•!> + f^)^^ where C, is the bulk viscosity, 77 is the shear 
viscosity, and the speed of sound is Vs = \JdP/de. 

^°The second counterterm, proportional to the boundary Ricci scalar, does not contribute to this order 
^'in AdS background, as opposed to the AdS-Schwarzschild case we consider, the on-shell action linear in fluctua- 
tions is zero. 
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To get the correct dimension for the stress tensor, we can restore the temperature dependence 

(r°)o = (27rT)^ x(-^)= -lir'NlT\ (118) 



E.2 Second order action 



Now expand (1581) to second order in the linearized on-shell fluctuations. To highlight the role of 
the second counterterm in (|58T ). and explain why one needs to go beyond the quadratic action given 
by [l26l . we first expand 

^{R-2A)+ [ ^^y{a + 2K) (119) 

M JdM 

and collect the {5gY terms. Since the fluctuations obey the linearized equations of motion, we 
are left again with a total derivative term (basically the first-order expansion of (|1 141 )). (The bulk 
term cancelled in the first-order expansion of the action (11191) because the background solved the 
Einstein equations. The bulk term cancels now because the fluctuations are on-shell.) After these 
considerations, this is what is left: 

1 
4 



u=0 



{d,g,,)5g;5g^^ - {d,g,,)g^' 5g^:,5g; - \{d^g^.)bg^'' bg^ 



-r\{d^g^.)g^'bglbgl \ + \abg^^bgl - abg^M""^ 



(120) 



Again, the boundary volume counterterm, with its a coefficient, removes a 1 /m^ divergence from 
the contact terms in (|120l) . besides contributing to the finite terms. Upon substituting the black 
hole background metric (|62] ). and 



rc' = -uVl - u^, J-gbdy = ^ ^ (121) 

V 



2 



in the second order expansion (11201 ) we arrive at 

1 f 1 



u=0 



U 



^\ {^liKf - + h% + - Ih^hi^ (122) 



where fi,^, p = 0, 1, 2, 3; i, j = 1,2,3 and Einstein summation convention was useco- 

However, while the action (II 191 ) is properly regularized and gives correct answers for the field 
theory stress tensor two-point functions up to terms linear in tu, k (and therefore suffices as long 



^^This agrees with (3.15) in |26|. An apparent discrepancy is resolved by noting that the derivative terms in [261 
include 95((/io)^), whereas we have d5{h'^hQ). This difference is reflected in a different coefficient of h'^hg. 
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as one is interested in linear hydrodynamic coefficients), it diverges as 1/m to quadratic order in 
uj,k. To remove this divergence a second counterterm is needed, proportional to the boundary 
Einstein-Hilbert action [|20l 1211 1221 1231 : 

-RAdS 



d-2 



fi'bdy-Rbdy (123) 
dM 



where we recall that for us = 4 and we set -RacIS = 2. It is easy to see that this second counterterm 
contributes to leading order only to terms quadratic in to, k. 



1 f VT 



u 



2 Ju=, 



d'bg^^ + d^dM', - d.dngi - d^dM^^ - b^AQ'K - d'djg;) Sg^ (124) 



since the curvature tensor is quadratic in derivatives. In (11241) indices are raised and lowered on 
the partial derivatives with the boundary metric %dy^i, = diag(— //u, 1/u, 1/u, l/u). Since there 
is one inverse metric per term in (I124] | we see that the contribution coming from this counterterm 
will be divergent as l/u. The conclusion is that this counterterm's job, to second order in fluctua- 
tions and second order in to, k, is only to remove divergences from (I122] |. without any finite term 
subtraction. At fourth order in A; we get a finite contribution, coming from terms linear in u in 
5gj. For completeness, we give below the leading order in to, k of (|124|) : 



1 

u=0 



(co' - k')[ihlf - hlhl] + ik'h', - hlco')ih\ + hi) 

-ukhiihi + hi) + ojk{hih\ + hlhl) + uj\h\)' + [hif] + e[{hif + {hi) 



(125) 

The second-order on-shell gravitational action S^'^^S is given by the sum of equations (11221) and 
(11251) multiplied by the gravitational prefactor N^/{2^7i''). 

E.3 Cubic action 

In writing the cubic action, it is helpful to start with the second-order Einstein-Hilbert plus cosmo- 
logical constant action 



5(^)5eh = S^''> [ ^{R - 2A) 
Jm 



M 

'v'\Ug^HgKL - \6g^Sgi) + iV^(5^;f 5/^^) 



V^Vm 

M 



-2Vj,{5g''^5g^) + dg^'^Vj.dg'J + 5/^^V^5^7^ 



(126) 



where the indices are being raised and lowered with the background metric (Jmn, and the deriva- 
tives are background covariant. In general, there is another contribution to the second-order action, 
which is proportional to the linearized equation of motion. Since 6gMN solves the linearized 
equation of motion, the only non- vanishing contribution is the total derivative given in (|126l) . In 
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expanding to terms cubic in the fluctuations, the Einstein-Hilbert plus cosmological constant action 
receives contributions from S^^^ (^'•^^^eh) as weU as from the expansion of the term proportional to 
the linearized equations of motion: 



MN 



(127) 



M 



A brute force calculation gives 



5(1) 5(2) c 



EH 



M 



M 



+ Ug'^V'\5g^p5g''^) + ^Jg^P^g'^^V'' 5g^^ 



l5g'^^5g^pVM5g^ 



ISg^Sg'^^'V^Sg^' - Ug^Sg^'^'V ^Sg^p + Ug'^^dgpQV ^Sg''^ 



MQ 



+6g''^6gpQVM6g''^ - 
-|5/^^V^(5/«5^pq; 



MP i 



"""^VNSgpQ - j-,6g^'^5gpQVN5g 



MN 



-,NP f^Q-r-jM 



(128) 



The Lmn tensor which appears in (11271 ) was introduced in [|30 



— \Lmn = S^^^Rmn — \5gMN{R — 2A) 

= -\{V^VK^gMN + ^M^N^gK - y^^M^gNK - N^gMx) + d6gMNR~KdS 

(129) 

where we recall that the background Ricci scalar is _R = 2{d + l)A/(c? — 1) = —d{d + l)-RAds ^^'^ 
(i = 4 for us. In terms of Lmn, the linearized equation of motion is Lmn — \gMNL = 0. 

In [f30l it was also shown that the second order variation of the equation of motion, |5(i) {Lmn — 
\gMNL) is expressed in terms of another tensor, Vmn'- 



- i5(i)L 



MN 



= V, 



MN 



5g^^{ypVQ5gMN + VmVjv^^pq - VpV 



M 



- ^p^n^Qmo) 



--(2Vp5(7^^ - V'^5g'p){VM6gNQ + ^N^gMQ - Vq5, 
+^[(Vm5^?''^)V^5(7pq + 2{V''5gl)Vp5gQN - 2(V^^ 



JMN) 

9m)^Q^9pn]- 
(130) 



Lastly, substituting into the second term of the cubic action (I127h . one finds [|30l 

|f5(i) / ^\{LMN-\gMNL)5g'''' = -\ f ^{yMN~\gMNV)5g''''.{\3\) 
Jm Jm 

The total derivative terms omitted from the cubic action in [|30l are those we give in (|1281 ). We do 
not throw away these terms, and keep their contribution (from both boundary and horizon) to the 
retarded momentum space 3-point stress tensor correlators. 



^^Lmn given in equation (2.6) in lf30ll has been further simpHfied by commuting two covariant derivatives, and 
subsequently using that the AdS curvature tensor is Rmpn^ = —{gMNSp — .9PA"^M)^AdS- Arutyunov and Frolov 
have further set RAdS = 1- 
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For the correlators we evaluate in Section 3.2, we only need cubic vertices of one tensor mode 
and two shear modes. This cuts down significantly the number of relevant terms from the cubic 
action (I127D . 

The Gibbons-Hawking term contribution to the cubic action is most easily evaluated by making 
explicit use of the gauge condition 5gM5 = and of the form of the unit normal vector to the 
boundary (11211) : 



The cubic action 5^^^S is given by the sum of ^^^^iSeh + '^'^^^'^gh and of the cubic expansion of the 
boundary counterterms, multiplied by the gravitational prefactor A^"^ / (2^7r^). 
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